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ABSTRACT 

We report observations of the CO J = 7 — > 6 transition toward the starburst 
nucleus of NGC 253. This is the highest-excitation CO measurement in this 
source to date, and allows an estimate of the molecular gas excitation conditions. 
Comparison of the CO line intensities with a large velocity gradient, escape 
probability model indicates that the bulk of the 2-5 xlO 7 M Q of molecular gas 
in the central 180 pc is highly excited. A model with T ~ 120 K, nn 2 ~ 4.5 x 
10 4 cm -3 is consistent with the observed CO intensities as well as the rotational 
H2 lines observed with ISO. 

The inferred mass of warm, dense molecular gas is 10-30 times the atomic gas 
mass as traced through its [C II] and [O I] line emission. This large mass ratio 
is inconsistent with photodissociation region models where the gas is heated by 
far-UV starlight. It is also not likely that the gas is heated by shocks in outflows 
or cloud-cloud collisions. We conclude that the best mechanism for heating the 
gas is cosmic rays, which provide a natural means of uniformly heating the full 
volume of molecular clouds. With the tremendous supernova rate in the nucleus 
of NGC 253, the CR heating rate is at least ~ 800 times greater than in the 
Galaxy, more than sufficient to match the cooling observed in the CO lines. 

Subject headings: ISM:molecules- galaxies:individual(NGC 253) — galaxies: ISM 
— galaxies:nuclei — galaxies: starburst — submillimeter 
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1. Introduction 

NGC 253 is a nearby (d ~ 2.5 Mpc; Mauersberger et al. 1996 ), nearly edge-on (i ps 78°, 
Pence 1981), isolated starburst galaxy with the starburst strongly concentrated at the nu- 
cleus. Based on the distribution of the 10-30 /im continuum, the bulk of NGC 253's starburst 
is confined to a ~ 60 pc region centered just south-west of the dynamical nucleus (Telesco, 
Dressel, & Wolstencroft 1993). The 1-300 fim luminosity within about 30" is 1.6 x 10 10 L 
(Telesco & Harper 1980) while the total IR luminosity of NGC 253 detected by IRAS is 
Lir w 2 x 10 10 L (Rice et al. 1988). There is a clear bar in the 2 MASS image of the galaxy, 
and the kinematics show evidence for orbits in a bar potential (Scoville et al. 1985; Das et al. 
2001, and references therein). The bar likely plays a role in transporting gas to the nucleus 
and thus supporting the starburst. Current estimates of the neutral gas mass in the central 
20" - 50" range from 2.5 x 10 7 M (Harrison, Henkel, & Russell 1999, hereafter HHR99) to 
4.8 x 10 8 M (Houghton et al. 1997). 

As is the case in many of the nearby starburst galaxies, much of the molecular gas 
in NGC 253 appears to be highly excited (Wild et al. 1992; Mao et al. 2000; Ward et al. 
2002). Observations of J = 4 -> 3 and J = 6 -» 5 transitions of CO (Israel & Baas 2002; 
Harris et al. 1991), as well as J = 1 — > and J = 3 — > 2 transitions of HCN (Paglione et al. 
1995, 1997) suggest that a large fraction of the molecular gas in the nucleus of NGC 253 has 
n Ha > 10 4 cm~ 3 and T ~ 100K or higher. Observations of higher-J CO transitions are critical 
to measure this warm molecular gas component, and probe its physical conditions. In our 
Galaxy, the warm, dense molecular gas conditions which excite the mid-J CO lines (J = 4-8) 
are typically associated with molecular cloud surfaces exposed to stellar far-UV radiation. In 
these photodissociation regions (PDRs), the UV dissociates and ionizes the gas as it enters 
the cloud and is absorbed by dust. The result is a layered geometry based on the extinction 
of UV; the progression from the exterior inward is ionized, neutral atomic, warm molecular, 
and finally cool molecular gas on the interior (Tielens & Hollenbach 1985; Wolfire, Tielens, 
& Hollenbach 1990; Kaufman et al. 1999). Several of these regions have been studied in 
star formation sites within the Galaxy (Harris et al. 1985; Schmid-Burke et al. 1989; Stacey 
et al. 1993), and they have been discussed in the context of starburst nuclei (Mao et al. 
2000; Carral et al. 1994). Other potential molecular gas heating sources are outflows and 
cloud-cloud collisions in which shocks convert dynamical energy to thermal energy of the gas, 
such as in the circumnuclear disk of the Galactic Center, and in the massive star formation 
sites W 51 and DR 21 (Jaffe, Harris, & Genzel 1987; Jaffe et al. 1989). Finally, it has been 
proposed that much of the molecular gas heating in starburst nuclei may be due to cosmic 
rays (Suchkov, Allen, & Heckman 1993, hereafter SAH). To better understand the energy 
balance of the molecular gas in the starburst nucleus of NGC 253, we have observed the CO 
J = 7^ 6 transition with our submillimeter imaging spectrometer, SPIFI on the JCMT. 
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2. Observations 



NGC 253 was observed in April 1999 during the commissioning run of the South Pole 
imaging Fabry-Perot interferometer (SPIFI) (Bradford et al. 2002) at the James Clerk 
Maxwell Telescope (JCMT) on Mauna Kea, Hawaii. The Fabry-Perot interferometer (FPI) 
was set to a velocity resolution of 70 km s _1 and was scanned over 620 km s" 1 with the 
center wavelength of the scan set to the redshifted (250 km s _1 ) CO J = 7 —> 6 rotational 
transition. The FPI scan length was measured with interference fringes of a He-Ne laser, 
and absolute frequency reference was an absorption spectrum of the J = 6 — > 7 CO transition 
from gas in a cell against a hot load. On this observing campaign the bolometer array was 
equipped with 12 functioning detectors arranged on a 4 x 4 grid with four blanks. Because 
of a figuring error in the polyethylene lenses used as a f-converter, the main beam as mea- 
sured with Mars had a size of 11.5" (FWHM), 1.6 times bigger than the diffraction-limited 
beam size of 7" per pixel. Mars was also used for flux calibration assuming a brightness 
temperature of 201 K (Wright 1976; Rudy et al. 1987). Flat-fielding of the array was accom- 
plished by rastering Mars around the array. For the period of the observations, the zenith 
atmospheric transmission was between 0.22 and 0.29, giving a transmission to the source 
between 0.05-0.12. This low and variable atmospheric transmission provided the greatest 
calibration uncertainty, which we estimate to be about 30%. For the observation of NGC 253 
the array was centered at the IR nucleus at R.A. = 0M5 m 06 s .0, Dec. = -25°33'40" (1950) 
((0, Opposition). We estimate our pointing accuracy to be 5". The data were obtained by 
chopping and nodding the telescope and the data presented here were obtained with 20 min- 
utes on-source integration. The 50 minute duration of the observations resulted in a field 
rotation of 18° about the center of the array, producing a total travel of 7" (±3.5" about the 
mean) for the corners of the array. Because this is comparable to the pointing uncertainty, 
and is substantial only for the corner pixels, we neglect the rotation and adopt the mean 
position angle. 



3. Results 

3.1. Comparison with Other Molecular Tracers 

Figure 1 shows the spectra obtained in this single array pointing toward the nucleus 
of NGC 253. The data are plotted in main beam brightness temperature units, correcting 
for beam efficiency, atmospheric transmission, and source coupling. Typical spectra peak at 
3-5 K, with a single spectrum showing a peak of 12 K; integrated intensities are between 
400 and 1200 K km s -1 . The integrated intensities in K km s" 1 are roughly half of the 
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CO J = 6 — > 5 intensities of Harris et al. (1991), when the J = 6 — > 5 antenna temperatures 
are corrected to a 10" main beam. The peak temperature of 12 K at v — 150 km s _1 in 
the spectrum just northeast of the nucleus has a similar peak above the Gaussian profile 
in the J = 6 — > 5 spectrum toward the nucleus. In addition to the bright spectra near the 
nucleus, we detect emission at 30-40% of the nuclear value at positions 10"-20" offset from 
the nucleus. Harris et al. also detect strong CO J = 6^5 emission at offset positions. They 
measure an integrated intensity 60% of the nuclear value as far as 14" offset, consistent with 
our observations. 

Interferometric observations of millimeter wave CS (Peng et al. 1996), CO (Canzian, 
Mundy, & Scoville 1988; Carlstrom et al. 1990), and SiO (Garcia-Burillo et al. 2000) transi- 
tions show evidence for a nuclear bar with major axis oriented 70° east of north. The velocity 
field of this bar shows a 100 km s^ 1 shift from the northeast to the southwest. Though the 
J = 7 — > 6 map does not have the resolution or sampling to spatially resolve the bar, it also 
shows the velocity shifting from 150-200 kms -1 in the northeast to 250-300 kms" 1 in the 
southwest, showing the overall kinematics of gas in the bar, in good agreement with the 
interferometric observations. The observed line widths are between 200 and 300kms _1 , con- 
sistent with the 150 to 240 kms -1 line widths observed in the lower- J CO lines with 12 to 
15" beams (Mauersberger et al. (1996); HHR99), convolved with our 70 kms" 1 Lorentzian 
instrumental profile. 



3.2. Model fit to the nuclear excitation conditions 

To make quantitative estimates of the physical parameters of the nuclear molecular gas 
we have compiled all CO rotational transitions observed toward the central 15" available in 
the literature to compare with predictions of an escape probability, large velocity gradient 
model for CO line excitation. In contrast with previous studies of the molecular material 
in starburst galaxies which invoke two components (HHR99) or complex PDR models (Mao 
et al. 2000), our modeling suggests that the observed 12 CO and 13 CO intensities can be 
reproduced with the bulk of the gas mass at fairly high excitation. The line ratios can be 
fit with a range in conditions from T = 100 K, nn 2 ~ 5 x 10 4 cm~ 3 to T > 300 K, nn 2 ~ 
2.0 x 10 4 cm~ 3 . Recent measurement of rotational H2 lines (Rigopoulou et al. 2002) suggests 
that the temperatures may be on the low end of the allowed range, and we thus we favor 
T ~ 100 K and n ~ 5 x 10 4 cm~ 3 . Combining the 13 CO transitions constrains the column 
density to be only slightly larger than that measured by HHR99 in a 23" beam. To reproduce 
the observed intensities with the LVG formalism requires large velocity gradients (dv/dz ~ 
10 — >• 80kms _1 pc~ 1 ), which are much higher than is traditionally used, but nevertheless 
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plausible for the starburst nucleus. Here we describe our modeling approach and the basis 
for these conclusions. 

Including our results, 12 CO has been observed toward the nucleus of NCG 253 in all 
transitions up to J = 7 — > 6 except for J = 5 — > 4, with a variety of beam sizes. Also, the 
optically thin isotopes 13 CO and C 18 have been observed in transitions up to J = 3 — > 2. 
Table 1 summarizes the 12 CO and 13 CO observations and the corrections to refer them to the 
15" beam. Measurements for all J = 3^2 and lower transitions are from HHR99. Coupling 
corrections, where necessary, are based on power laws. The J = 2 — > 1 lines have been observed 
with both 12" and 23" beams, and the resulting power law is I oc Q~ - 76 . For J = l— ►(), there 
are no single dish observations with a beam smaller than 23", but based on the centrally 
peaked J = 3^2 / J = 2^1 ratio map, the J = 1^0 emission is probably more extended than 
J = 2— >1. We adopt a power law of I oc 9~ 0A , but estimate a 20% uncertainty in the derived 
intensity, which is sufficient to accommodate a power law exponent between and 0.76. For 
J = 3^ 2, the intensities are measured from the published spectra in 15" beams. For J = A—> 3 
there are published measurements with both a 10.6" beam (Israel et al. 1995) and a 22" 
beam (Israel & Baas 2002). Allowing these two measurements to define a power law gives 
I oc 6~ im , and we use this to determine the intensity within 15". The J = 6^5 measurement 
is from Harris et al. (1991) in T^ units, and is corrected to 15" according to their stated 
8"/30" composite beam profile. The J = 7 — > 6 observation from this work is referred to 
an 11.5" main beam based on observations of Mars, with a small correction for the 15.8" 
planetary disk. The estimated uncertainties are also included in Table 1. The beam filling 
corrections become the dominant source of uncertainty in the J = 1 — > (described above) 
and J = 2 — > 1 lines, we estimate this to be about 15%. For the J = 3 — > 2 transition in a 15" 
beam, we estimate 20% for the overall uncertainty including systematics and statistics. For 
the J = 4 — > 3 and higher submillimeter transitions, the uncertainties in the measurements 
themselves are substantial, as much as 30%. 

The LVG model calculates intensities of rotational CO lines given a molecular hydrogen 
density, a temperature, and a velocity gradient. The calculation is similar to that used in 
other multi-line analyses, iteratively solving for the level populations up to J = 15. We 
use the escape probability j3 — (1 — e~ T ) jr of Goldreich, & Kwan (1974) appropriate for 
spherical clouds undergoing gravitational collapse. H 2 - CO collisional excitation rates are 
those of Flower (2001), increased by 21% to account for collisions with He (McKee et al. 
1982; Viscuso & Chernoff 1988; Flower & Launay 1985; Flower 2001; Schinke et al. 1985). 
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Table 1. CO observations toward the NGC 253 nucleus 



Transition 


v [GHz] 


Obs. Beamsize 


/ T MB dv 


Corr'n to 15" 


Int. in 15" 1 


± % 


12 CO J = 1^0 


115.3 


23" 


920 


1.2 


1.72 x 10~ 6 


20 


12 CO J = 2^1 


230.5 


12" 


1060 


0.84 


1.37 x 10" 5 


15 


12 CO J = 3^2 


345.8 


15" 


998 


1 


6.5 x 10" 5 


20 


12 CO J = 4^3 


461.4 


22", 10.6" 


1019, 2150 


1.48, 0.70 


1.50 x 10" 4 


30 


i2 C J = 6^5 


691.5 


8" / 30" 


861 


1.59 


4.6 x 10- 4 


30 


i2 C J = 7^6 


806.7 


11.5" / 60" 


1370 


0.677 


5.0 x 10~ 4 


30 


13CO J = 1^0 


110.2 


23" 


80 


1.2 


1.32 x 10~ 7 


20 


13 CO J = 2^1 


220.4 


12" 


134 


0.84 


1.24 x 10~ 6 


15 


13 CO J = 3^2 


330.6 


15" 


145 


1 


5.4 x 10~ 6 


20 



Note. - - 1 Integrated intensity in ergs -1 cm -2 sr" 1 . The intensities for J = 3^ 2 and lower 
are taken from Harrison, Henkel, & Russell (1999). J = 4^3 is from Israel et al. (1995); Israel 
& Baas (2002), J = 6^5 from Harris et al. (1991), and J = 7^6 from this work. The couplings 
and uncertainties are described in the text. 
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3.2.1. 12 CO Intensity Ratios 

A sensible approach for comparing an LVG model with observations is to use the line 
intensity ratios to determine the physical parameters (temperature, density and opacity) 
while the absolute line intensities determine the total column density. In addition to the 
12 CO intensity ratios, the optically-thin isotopes can be useful in constraining the 12 CO 
optical depths, which are greater than unity for the lines of interest. In order to concentrate 
on the nuclear gas, we do not put much emphasis on the J = 1 — ► transition. Since it 
connects the ground state and and the J = 1 level is only 6 K above ground, the J = 1 — > 
transition traces any molecular material in the beam. A related and more confusing problem 
is the potential for self-absorption by quiescent material along the line of sight but exterior 
to the nucleus, a particular concern for nearly edge-on systems such as NGC 253. This may 
be the reason that any model which fits the J = 2 — > 1 and higher transitions over predict 
J = 1 — ► 0. In an approach similar to that of HHR99, we invoke a small diffuse component 
which is negligible energetically, but which emits primarily in J = 1 — > to account for its 
observed intensity. 

For a given LVG calculation, the parameters which determine the line intensity ratios 
are the temperature, density, and column density of CO per unit velocity interval, Nco/dv, 
called the opacity parameter. Nco/dv is a combination of density, velocity gradient of the 
gas, and CO abundance, N co /dv = X C on H2 (dv/dz) _1 where dv/dz is the velocity gradient 
in km/spc -1 . We have assumed a CO abundance Xco of 8 x 10~ 5 , the value measured 
for galactic clouds (Freking, Langer, & Wilson 1982). While the metallicity in the nucleus 
of NGC 253 is as much as 3 x solar (Mauersberger et al. 1996), the large observed [CI] 
intensities suggest that, relative to the Galaxy, the gas-phase carbon budget includes a 
much larger neutral carbon fraction, and smaller CO fraction (Israel & Baas 2002), thus the 
CO abundance is not clearly much higher in spite of the large metallicity. Intensity ratio 
calculations for a higher CO abundance are identical to those with the canonical Xco value, 
but with the velocity gradient decreased by the same ratio. Raising the CO abundance or 
lowering the velocity gradient increases the the optical depth in the lines. Figure 2 plots the 
results of LVG calculations of diagnostic line ratios as a function of gas temperature, density, 
and velocity gradient. We have found that the J = 7^ 6 to J = 3^ 2 intensity ratio is a good 
measure of the overall excitation, and that matching this ratio typically satisfies all the lines 
above J = 1 — > to within the observational uncertainties. 

As can be seen from the solid contours in Figure 2, the 12 CO J = 7 ^ 6 to J = 3 ^ 2 
intensity ratio of 7.7 means that the gas is of high excitation. Because the run of line 
intensity with J is increasing up to J = 7 —> 6, there is no unique solution, nor even upper 
limits to the temperature or density. Lower limits, however, are strong. For temperatures less 
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than 300 K, the density must be greater than 2 — 3 x 10 4 cm~ 3 . A temperature less than 150 
K requires a density greater than ~4x 10 4 cm~ 3 , clearly different from the typical Galactic 
molecular cloud. To first order, these results are independent of the velocity gradient and 
CO abundance because the lines are optically thick, even at velocity gradients as large as 80 
kms -1 pc _1 . 



3.2.2. 13 CO Intensities 

If their abundances relative to the main isotope are known, the intensities in the 13 CO 
transitions can be used to determine the 12 CO optical depths. This is the approach taken by 
HHR99 in their study of the nuclear regions of NGC 253 with the lower-J lines. The key to 
their analysis is the assumption of a 12 CO to 13 CO abundance ratio of 40-50, originating from 
analysis by Henkel et al. (1993); Henkel & Mauersberger (1993). Our results are consistent 
with a ratio on the upper end of this range, and we have chosen 50 as a value to use in 
the model plots. Included in Figure 2 as dashed lines are the 12 CO to 13 CO J = 3 — > 2 
intensity ratios predicted by the LVG models. The observed ratio of 12 is plotted as a 
thick dashed line. To correctly obtain the line intensity ratios between the two isotopes 
requires running the full LVG model for the optically thin species, since radiative trapping is 
important in determining the level populations of 12 CO but is much less so for the optically 
thin species. Because the 13 CO transitions are optically thin for all relevant conditions, 
the values plotted in Figure 2 can be shifted to accommodate a different isotope ratio. In 
the density-temperature plane, a larger isotope ratio is considered by choosing a line ratio 
contour with a value decreased from the observed 12 by the same factor. 

In the panels of Figure 2, the temperature and density at which the two observed line 
ratios meet is the regime of possible solutions. It is evident that small velocity gradients are 
unable to generate the observed line ratios. At these small velocity gradients, the conditions 
which excite upper-J lines generate large optical depths in the 12 CO lines, and their fluxes 
are too small compared with the optically thin isotope. Increasing the velocity gradient 
does provide possible solutions, again, with high excitation conditions in both temperature 
(T>100 K) and density (nn 2 > 4 x 10 5 ). Though unconventional in LVG analyses, the large 
velocity gradient is plausible for the starburst nucleus, and is discussed in Section 3.2.4. 
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3.2.3. CO Column Density 

Even though the 12 C0 lines are optically thick, their observed intensities put lower 
limits on the CO column density independent of the optically thin isotopes. The intensity 
in a line is determined by the product of column density in the upper level and the optical 
depth factor 1 — e~ T , hence a lower optical depth (larger velocity gradient) generates more 
emission for a given CO column density. The curves in Figure 3 show the column density 
required by the LVG model to reproduce the observed intensities. For each velocity gradient, 
the curve is parametrized in terms of temperature, with the density determined according to 
the observed 12 CO line ratio shown in Figure 2. Examination of Figure 3 indicates that the 
column density of warm gas must be substantial. For example, if the gas has a temperature 
below 200 K and velocity gradient below 40 kms -1 pc -1 , the column density must be larger 
than 3 x 10 18 cm -2 , a lower limit which is comparable to the 3.5 ±1.5 x 10 18 cm -2 estimated 
in the central 23" with the J = 3^ 2 and lower transitions (HHR99). 

To improve upon the lower limits provided by the 12 CO transitions alone, we incorporate 
the 13 CO intensities. The boxes in Figure 3 show the column densities required by the LVG 
model under the conditions corresponding to the intersections of the two diagnostic line 
ratio curves in the panels of Figure 2. Boxes are plotted for values of dv/dz = 10, 20, 40 
and 80 kms _1 pc _1 , as velocity gradients less than 10 kms _1 pc _1 do not allow reasonable 
solutions. The sizes of the boxes correspond to a ±30% uncertainty in temperature, and 
an corresponding uncertainty in column density given by local slope of the T-Nco relation. 
The locus of these regions indicates that the CO column density is between 3-7xl0 18 cm -2 , 
for all allowed values of the excitation. This is consistent with the 12 CO upper limits, and 
is comparable to the value derived by HHR99 for the central 23". 

The results of our LVG analysis can be summarized as follows: the molecular gas in 
the central 180 pc of NGC 253 is much more highly excited than that traced with only the 
lower-J lines in a larger 280 pc beam. The gas is warmer than 100 K, and has a density 
greater than 2 x 10 4 cm -3 . The 12 CO lines have optical depths not larger than 2-4, requiring 
velocity gradients > 10 kms -1 pc -1 , larger than traditionally employed in LVG studies. Even 
larger velocity gradients (> 20 kms _1 pc _1 ) are required if the temperature is constrained 
to be below 200 K. The CO column density is well-determined by the 13 CO lines: N C o =3- 
7xl0 18 cm -2 , only slightly larger than the values estimated in a larger beam. Finally, while 
the 13 CO J = 1^0 transition requires an additional component of cool and/or diffuse gas, 
this component is small by mass, negligible energetically, and may not even be associated 
with the nuclear region. There is no reason to invoke a large component of cool gas as would 
be required, for example, in a study of the central regions of the Milky Way. 
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3.2.4- The Large Velocity Gradient 

Previous LVG studies of NGC 253 and M82 have employed a velocity gradient of around 
1 kms _1 pc _1 , which is approximately the ratio of the total observed velocity width to the 
size of the observed region (HHR99; Mao et al. (2000)), and comparable to the values used 
in studies of Galactic molecular clouds. This is a lower limit to the true velocity gradient in 
the material, appropriate if it were only subject to the central potential of the galaxy with 
no local motion. In reality, clouds and clumps will be subject to their own gravity as well 
as that of the stars around them. The escape velocity of a self gravitating cloud of radius 
r, and density p determines a velocity gradient given by dv/dr ~ a/87t/3Gp, which gives 10 
kms _1 pc _1 , for a density of nn 2 ~ 4.5 x 10 4 cm -3 . Models on the low end of the allowed 
temperature range, e.g. T ~ 120 K require a velocity gradient of order 40 kms -1 pc -1 , a 
factor of 4 larger than this escape velocity. This is not implausible; if the star-forming clouds 
in NGC 253 are similar to the Orion molecular cloud core, then one might expect the stellar 
mass density to be about 4 times the gas density (O'Dell 2001), which would double the 
effective escape velocity. Furthermore, Carral et al. (1994) find the pressure of the ionized 
gas in the central 45" of NGC 253 is ~ 7 x 10 6 K cm~ 3 - larger than the thermal pressure 
within our modeled molecular gas. This high pressure ionized gas has a large volume filling 
factor (0.01 to 0.3) within their much larger (45") beams. It is thus quite plausible that 
within our 15" modeled region the ionized gas place a important role in cloud confinement, 
further increasing the velocity dispersion in the gas. In short, the large velocity gradients 
required for a good model fit are very plausible given the dynamics in the nucleus of NGC 253. 



3.3. Comparison with rotational H 2 lines 

Because the J = 7—> 6 transition is not high enough to see the turn-over in CO excitation, 
the CO analysis does not provide a good upper limit to the gas temperature. Fortunately, we 
can use the ISO H2 pure rotational line observations to provide an important temperature 
constraint. Rigopoulou et al. (2002) have recently reported ISO SWS observations of the pure 
rotational lines of H 2 from the central 14" x 20" regions of NGC 253. These quadrupole H 2 
lines have low radiative rates, so they are easily thermalized and optically thin. In principle 
the lines can be used to measure temperature and column density. However, even the lowest 
H2 rotational transition energy is 515 K, thus the observed line emission is heavily biased 
toward the very warmest regions (which may be negligible by mass), and each observed 
transition is an exponential function of gas temperature. Not surprisingly, the line intensities 
are not consistent with a single excitation temperature. Based on the S(0), S(l), S(2), 
S(5), and S(7) lines, Rigopoulou et al. (2002) conclude that there must be a substantial 
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mass of warm, dense gas in the nucleus of NGC 253, broadly consistent with our results. 
More specifically, the modeled molecular hydrogen column density of N Ha = N co /X co ~4- 
8xl0 22 cm~ 2 is consistent with the S(0) data presented by Rigopoulou et al. (2002) for 
gas temperatures in the range of 95-120 K. We therefore favor models on the low end of 
the temperature range allowed by CO analysis. Because of the exponential dependence on 
temperature, this H2 temperature constraint is a rather weak function of H2 column density, 
hence also a weak function of CO column density and CO abundance. 

We combine this H2 result with the parameter ranges allowed by the LVG analysis, and 
the physical constraint of a small velocity gradient to arrive at our preferred model, shown in 
Figure 4. This model uses dv/dz = 40 kms -1 pc _1 , toward the bottom of the range allowed 
by the low temperature and the observed 13 CO intensity. As described above, the lower 
velocity gradients are more consistent with gravitationally bounds clouds. The modeled 
conditions are T = 120 K, n H2 = 4.5 x 10 5 cm~ 3 , and the CO column density required to 
generate the observed line intensities is 5.3 x 10 18 cm~ 2 . At this temperature, about 40% 
of the observed H2 S(l) flux is generated in the bulk component, and the contributions to 
the higher-J H2 lines (which require excitation of levels with E > 1500 K) are less than 
10% of the observed line fluxes. As an estimate of the uncertainties, we note that even with 
the temperature constrained, a larger velocity gradient still allows a smaller column density. 
Decreasing the velocity gradient, or increasing the CO abundance can permit a larger column 
density, but this quickly becomes inconsistent with the 13 CO observations unless the 13 CO 
relative abundance is decreased below our adopted 1/50. Combining all the uncertainties 
in the abundance ratio, modeling, and observations, we estimate an uncertainty of 40% for 
the final adoped column density: Nco ~ 3.5 — 7 x 10 18 cm~ 2 . Our model parameters are 
summarized in Table 2. 



3.4. Warm Gas Mass 

With the assumed CO abundance relative to H 2 of 8 x 10~ 5 , the CO column density 
corresponds to a mass of excited molecular gas in the central 15" of 2.5-5. lxlO 7 M . In 
reality, the CO abundance is not well-known. The abundance may be larger than the canon- 
ical Galactic value, given the 3x solar metallicity (but see discussion above), and this would 
result in a smaller gas mass by as much as a factor of 3. Based on an assumed carbon gas to 
dust ratio and the measured fractions of gaseous carbon bearing species, HHR99 infer a CO 
abundance relative to H 2 of 2.2 x 10~ 4 , that is, about the factor of 3 larger than our adopted 
value. With this abundance, the gas mass would be smaller: 0.9-1.8 xlO 7 M . We do note 
a larger CO abundance demands even larger velocity gradients (more than 30 kms _1 pc _1 ) 
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to allow working solutions. For comparison, the most recent estimate of the gas mass based 
on 1.3 mm dust continuum is 4.4 x 10 7 M o (HHR99, Kriigel et al. (1990)). Since the 1.3 
mm continuum has a half-power size of only 16" x 11" much of this flux arises within our 
within our 15" beam, and traces the same material. The agreement with our derived gas 
mass suggests that while our CO abundance may be slightly high, it is not off by more than 
a factor of 2. 



3.5. Suitability of Multiple-Component Models 

It is traditionally thought that mid-J CO emission originates in the PDR surfaces of 
molecular clouds exposed to UV starlight. In this scenario, the warm excited gas is separate 
from the bulk of the molecular gas, and usually represents a small fraction of the total mass. 
Observations of several PDR regions in the Galaxy have motivated modeling which includes 
detailed PDR chemistry and radiative transfer to predict the line emission. Koster et al. 
(1994) have constructed detailed models of PDRs on the surfaces of small clumps to predict 
12 CO and 13 CO line intensities in star-formation regions. Relative to a simple plane-parallel 
model, the clumpy geometry allows an increase in local gas density, and allows UV photons 
to penetrate to a greater fraction of the gas, increasing the overall molecular excitation level. 
While such models are successful in reproducing the intensity ratios between 12 CO lines, 
they fail to account for both the large integrated intensities and the large 12 CO to 13 CO 
intensity ratios observed in a starburst nucleus such as NGC 253. For example, Koster et al. 
(1994) model slabs with thickness given by Ay = 10, illuminated from both sides, including 
chemistry of both 12 C- and 13 C-bearing species. The ratios among the 12 CO lines are well- 
fit with n H2 = 10 5 cm~ 3 , xuv = 10 5 Habing (= 1.6 x 10 2 ergs _1 cm -2 ), but to generate 
the observed J = 7 — > 6 flux in the central 15" of NGC 253 with a superposition of these 
regions would require between 5 and 20 times the CO column density estimated with our 
LVG analysis. The 13 CO lines show the same effect - the predicted J = 3 ^2 12 CO to 13 CO 
intensity ratio is around 3, compared with the observed value in 15" of ~ 12 (HHR99). 
Thus the existing PDR models are not appropriate for the global properties of a starburst 
nucleus. The key problem is the high optical depth in the PDR models, resulting from the 
small assumed velocity width per gas column. The effective velocity gradient of the clumpy 
PDR model is around 9 kms -1 pc _1 , too low as we have discovered in Section 3.2. 

To consider a PDR scenario with the appropriate optical depths in the lines for NGC 253, 
we have examined LVG models with two distinct gas components. Though this approach is 
not as rigorous as a full treatment of the chemistry and UV penetration used by the Galactic 
PDR models, to first order it accounts for the combination of warm UV-exposed surfaces and 
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cool cloud interiors which together would generate the observed CO emission. For NGC 253, 
we have found that while it is possible to fit the observed intensities without the full mass of 
gas being highly excited, the warm fraction must be substantial. To reproduce the observed 
J = 6^ 5 and J = 7^ 6 transitions, at least 30% of the total column density estimated with 
the 13 CO ratio must have high excitation conditions. Furthermore, the curves in Figure 3 
shows that this scenario requires an extremely large velocity gradient (~100 kms _1 pc _1 ) 
and/or very high temperatures (~300 K) in this 1/3 warm fraction in order to reproduce 
the observed intensities. Moreover, as we outline in Section 4.1.1, even with only 1/3 of the 
total column density accounting for the high-J line emission, the amount of warm molecular 
gas would still be 7-10 x greater than the PDR-associated atomic gas, so the UV heating is 
not a viable origin of the observed CO emission. 



4. Discussion 

4.1. Energetics of the Molecular Gas in NGC 253 

Independent of the details of the excitation and radiative transfer model, the general 
shape of the CO intensity vs. J curve measures the luminosity of the molecular gas in the 
central 15" of NCG 253. Taking our favored model plotted in Figure 4 as a guide, the total 
flux in the CO lines from the modeled region is 1.5 x 10 -11 ergs s _1 cm -2 . This is 25% of 
the [C II] 158 fim and [O I] 63 jum intensities from the entire inner 45" (Carral et al. 1994), 
and a substantial 1.4 x 10~ 4 fraction of the total IR luminosity of the galaxy. At a distance 
of 2.5 Mpc, this molecular gas flux corresponds to a luminosity of 2.9 x 1O 6 L , giving a 
luminosity to mass ratio of ~ 0.1 L©/M for the gas containing CO based on the mass 
estimates of Section 3.4. The large luminosity in the CO lines is a simple consequence of 
the fact that such a large fraction of the gas mass is highly excited. In what follows, we 
examine potential heating sources to generate the large luminosity. Simple arguments show 
that traditional heating sources for warm molecular gas - far-UV photons and cloud-cloud 
collisions - are unlikely to contribute substantially. In contrast, the enhanced cosmic ray flux 
of the starburst nucleus can easily match the observed cooling in the CO lines, and is likely 
responsible for the high excitation conditions inferred for the bulk of the molecular ISM. 

For completeness we note that the CO lines are a good measure of the total cooling of 
the molecular gas. The five H 2 lines observed by Rigopoulou et al. (2002) suggest that all 
the rotational lines combine to produce a flux of at most 1 x 10~ n ergs s _1 cm~ 2 , less than 
the CO lines. Another potential coolant is C°, which can have an abundance comparable 
to that of CO. However, it only has two transitions each of which could be comparable to 
a single CO line if optically thick and thermalized; thus C° cannot be a dominant coolant. 
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Also, we note that though the modeling described above carries uncertainties as great as 
40% in the derived parameters, the discussion below typically concerns factors of 3-30, and 
the general arguments hold independent of the modeling details. 

4-1.1. Stellar UV Heating 

In Galactic starformation regions, warm CO is typically produced by UV impinging on 
the surfaces of molecular clouds. As it penetrates a cloud, the UV is attenuated, producing 
in turn ionized, atomic, and warm molecular layers before it is extinguished (Tielens & 
Hollenbach 1985). The thermal pressure implied by our model fit to the CO lines is 5.4 x 
10 6 Kcm -3 , similar to the pressure derived with atomic PDR line ratios of 2.9 x 10 6 Kcm~ 3 
(Carral et al. 1994). This is consistent with a model in which the PDR gas forms the surface 
of the warm molecular clouds. The penetration of UV into a cloud is determined by the dust 
opacity, proportional to the gas column density according to 

A v = 2Ah 2 /1.9 x 10 21 cnT 2 . (1) 

Here we have assumed solar dust abundance. (This is not unreasonable — though the metal- 
licity in the nucleus of NGC 253 may be 3x solar (Mauersberger et al. 1996), there is 
evidence for substantial destruction of dust grains in an unusually large gas-phase Si abun- 
dance (Carral et al. 1994)). In standard Galactic PDRs, the first ~ 3 magnitudes from the 
surface of the cloud are photodissociated and cool through the atomic lines, while the warm 
molecular gas exists from A v ~ 3-6. To first order, the masses of these two components are 
comparable. Yet in the nucleus of NGC 253, Carral et al. (1994) measure 2.4 x 1O 6 M of 
photodissociated gas traced by the 158 /im [CII] line in a large (45") beam. Scaling by the 
observed molecular and far-IR continuum morphology (Canzian, Mundy, & Scoville 1988), 
(Telesco, Dressel, & Wolstencroft 1993), we estimate only 56% of the photodissociated gas 
traced with the large beam falls within our 15" beam. Therefore, the photodissociated gas 
mass is about a factor of 20-30 smaller than the 2-5 x 10 7 M of warm molecular gas. There 
is no reasonable geometry in which two layers of comparable thickness can have volumes 
different an order of magnitude. It is possible that the distribution of heating between the 
molecular and atomic components could be made to favor the molecular gas with soft or 
weak UV fields relative to standard PDR templates. However, the conditions derived with 
the atomic lines are consistent with standard photodissociation region (PDR) models - the 
UV field can be understood as arising from O 7.5 stars with an effective temperature of 
34,500 K (Carral et al. 1994; Engelbracht et al. 1998), so scenarios invoking soft or weak UV 
fields are unlikely. 

As pointed out in Section 3.5, as little as 30% of the total CO could provide the line 
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Table 2. Dense Molecular Gas Properties in the Central 180 pc of NGC 253 



Parameter Value 

Adopted 
CO abundance (X C o) 
12 CO / 13 CO abundance ratio 

Modeled 

NcO, beam, dense 

Gas density (nn 2 ) 
Temperature (T) 
Velocity gradient 
Dense gas mass 

Derived with best model 
Dense gas filling factor ((f>v) 2.6 x 10 

Cloud radius 2.0 pc 

Number of clouds 330 



8 x 10" 5 
50 

3.5-7xl0 18 cm" 2 
3-6 xlO 4 cm" 3 
100-140 K 
> 30kms _1 pc -1 
2.5-5 xl0 7 M m 



-16- 



luminosities in a warm component, emitting more efficiently lower opacity in the lines. The 
remaining 70% of the molecular gas would then be in a cool, more diffuse component, possibly 
the UV-shielded cloud cores. This scenario is attractive because UV-shielded regions are 
required to host the HCN observed toward the nucleus (Paglione et al. 1995), which has a 
low dissociation energy (~ 5 eV). However, a factor of 3 reduction in warm CO-containing 
gas mass would still leave a factor of 7-10 more warm molecular gas than atomic gas. Thus 
standard PDR models which invoke stellar far-UV photons to heat the surface of molecular 
clouds are not a good match for the heating of the bulk of the molecular gas in the nucleus 
of NGC 253. 



4-1.2. Dynamical Heating 

Another potential explanation for the large molecular luminosity is dynamical energy 
sources. By heating molecular gas over large velocity widths, shocks can dramatically in- 
crease the efficiency of CO radiation. One source of shocks are outflows from young stars 
or star associations, examples in our Galaxy include the W 51, W 49 and DR 21 regions 
(Jaffe, Harris, & Genzel 1987; Jaffe et al. 1989), in which bright CO J = 7 —> 6 emission is 
observed with ~20-50 kms -1 linewidths. In these sources, the CO luminosity in the broad- 
line shock-heated outflow components dominates that of the narrow-line PDR component. 
Another clue is the bright SiO observed in the central 150 pc of NGC 253 (Garcia-Burillo et 
al. 2000); the large implied abundance of SiO is attributed to outflows from massive young 
stars in a manner similar to the Galactic regions W 51 and DR 21. However, the comparison 
cannot apply globally in the nucleus of NGC 253. First, the excitation conditions in these 
Galactic shock regions are much higher that what we observe for the average properties of 
the nucleus of NGC 253. Observations of CO J = 2^1, J = 7^6 and J = 16^15 toward W 51 
IRS 2, for example, indicate that the gas temperature are 500-1000 K, unlikely for the bulk 
of the gas in NGC 253 given the rotational H2 observations suggesting lower temperatures. 
More importantly, in NGC 253, the near-IR H2 emission is characteristic of UV fluorescence, 
rather than the thermal spectrum produced by shock heating (Engelbracht et al. 1998) found 
in the Galactic outflow sources. 

Shocks can also be generated with large-scale gas motions: supernova blast waves or 
spiral density waves. The nuclear starburst in NGC 253 is approximately bounded by an 
inner molecular ring with diameter 120 pc which corresponds to the inner Inner Linblad 
Resonance (ilLR) of the galaxy's bar. The large non-circular motions observed (as much as 
40-50 kms -1 ) are attributed to gas accumulation and orbit crowding (see Buta & Combes 
(1996)), in the deep potential of the stellar mass, and it is plausible imagine that shocks from 
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cloud collisions would be common. However, the difficulty again is that if the majority of the 
CO luminosity is to come from cloud collision shocks, they must be slow (v ~ 10 km s _1 ), 
non-dissociative shocks in order to reproduce the observed CO line ratios, and not excite the 
ro-vibrational H 2 lines with a thermal spectrum. At these velocities, the collision timescale 
and energy available for collisions between clumps becomes small. Considering the molecular 
gas to come from assumed r~ 2 pc, nn 2 = 4.5 x 10 4 clumps, the mean separation of clumps 
A is around 25 pc. With a relative velocity of 10 kms -1 collisions will occur on timescales 
of 

r ~ A 3 /(^v) ~ 1 x 10 8 yr (2) 

with an available energy E = |M c i ump t> 2 . Such collisions, then, could in principle result in a 
luminosity to mass ratio of 

Mj dyn ™^~ 8Xl ° M~e (3) 

which is much smaller than the mass to luminosity relation implied by the model fit to 
the CO intensities (~ 0.1). Lower densities and smaller clumps sizes mean more numerous 
clumps and a smaller mean separation. This can increase the luminosity by a factor of 
only 10-20, but it is still several orders of magnitude too small. Thus it is unlikely that a 
substantial fraction of the CO luminosity in the central ~100 pc may be due to shocks from 
cloud-cloud collisions. Again, while higher-velocity shocks could provide much more energy, 
and the ~40 kms -1 non-circular velocities are suggestive that such shocks could exist, the 
H2 ro-vibrational lines are inconsistent with any significant heating in fast shocks. Thus 
while some amount of shock heating is required to produce the high SiO abundance, it is 
not likely to be a dominant factor in heating the bulk of the molecular gas. 



4-1.3. Cosmic Ray Heating 

It is widely believed that low energy cosmic-ray protons dominate the heating of the 
quiescent, UV-shielded molecular cloud cores in the disk of the Galaxy, keeping the gas 
temperature above that of the CMB (Goldsmith, & Langer 1978, and references therein). A 
much larger cosmic ray density generated by massive star formation has also been invoked 
to explain the 50-150 K measured molecular gas temperatures in the central regions of M82 
(SAH). Since the H 2 cosmic ray dissociation cross sections are so small (~3x 10 _26 cm 2 ), 
cosmic rays easily penetrate molecular cloud cores. The heating rate for cosmic rays can 
be characterized by a ionization rate ( p per H 2 molecule, and an average AQ deposited as 
thermal energy per ionization, that is: 

Xn 2 = ( P AQ. (4) 
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Goldsmith, & Langer (1978) adopt values for AQ between 17 eV and 20 eV. For the ionization 
rate we solve Eq. 3 of SAH which scales the cosmic ray density from that of the Galactic 
disk according to the supernova rate per unit area tp and a typical wind velocity v by which 
cosmic rays are removed from the disk of a galaxy: 

Cp,N253 n cr N253 '0N253 "^Galdisk /_>. 

-z = oc . (5) 

CpGal n cr> Gal WGal V W) N253 

The Galactic ionization rate C P ,Gai is between 2 x 10 _17 s _1 (Goldsmith, & Langer 1978) and 
7 x KT^s- 1 (van Dishoeck, & Black 1986), (see also Doty et al. (2002)). For the supernova 
rate in NGC 253 ^253, we take the lower limit derived for the central 360x90 pc of 0.1 
yr _1 (Antonucci & Ulvestad 1988), giving 3yr _1 kpc~ 2 . To compare with the Galaxy, we 
take a rate of 5.8 per century over the entire disk (van den Bergh & Tammann 1991, and 
references therein), or 8 x 10~ 5 yr _1 kpc~ 2 . The diffusion velocity in the Galactic disk is 
lOkrns -1 (SAH), and the superwind velocity in NGC 253 is ~500 kins" 1 , (Strickland et al. 
2002). The ionization rate can then be scaled from that of the Galaxy: 

Cp,N253 = 750 Cp, Gal = 1.5 - 5.3 x lO^V 1 . (6) 

With a AQ of 20 eV per ionization, the power deposited per H 2 molecule is then 5- 
18 xl0~ 25 ergs _1 . This estimate represents a lower limit to the cosmic ray heating, yet it is 
very close to the measured luminosity per H 2 molecule in all of the CO lines, 7 x 10~ 25 ergs -1 
with our assumed CO abundance. 

To check that the increased CR flux does not destroy H 2 and other molecular con- 
stituents, we compare our derived ionization rate with Eq. 10 of SAH. This requires that 
the ionization rate be less than the H 2 formation rate: 

C P N253 <1.5xl0- 18 nT 1 / 2 , (7) 

where n is the total hydrogen density. For our parameters, this requires CpN253 < 7.4 x 
10 _13 s _1 , a factor of 50 from our lower limit of C P N253- 

It is conceivable the increased CR density would dissociate the CO in the cloud core. 
Farquhar, Millar, & Herbst (1994) have constructed a chemical model of dark clouds to 
investigate the effects of cosmic rays ionization on the ion-neutral chemistry. For steady 
state models, with cloud density ~ 2 x 10 4 cm~ 3 and visual extinction A v ~ 30 they find CO 
remains the dominant reservoir for carbon for the entire range of cosmic ray ionization rates 
investigated, ( = 1.3 x 10~ 19 -^8x 10 _14 s _1 . The neutral carbon abundance grows sharply 
as C grows to 2 x 10~ 14 approaching 25% of the CO abundance at the largest cosmic ray fluxes 
investigated. These results are consistent with the large C°/CO abundance ratios observed 
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in the nuclei of both NGC 253 (Harrison et al. 1995), and M82 (Schilke et al. 1993). Since 
cosmic ray heating provides the necessary molecular gas luminosity and naturally accounts 
for the large observed neutral carbon abundance, while other heating sources are insufficient 
by large factors, we favor favor the cosmic rays as the dominant heating source for the bulk 
of the molecular gas in the nucleus of NGC 253. Our analysis is fairly general, and we note 
that similar arguments could be made for other starburst nuclei where the rate of supernovae 
explosions approaches 0.1 per year in a ~100 pc region. 

4.2. Cloud Size, Relation to Atomic Gas 

The total gas mass in a region compared with the derived density determines a volume 
filling factor of clouds. For our modeled dense component, the result is 0y = 2.6 x 10~ 3 . The 
physical size of the clouds, however, is not constrained by our observations. The existence of 
HCN in dense cores (Paglione et al. 1995) indicates that many clouds must be large enough 
to shield their interiors from UV, approximately A v > 10 for the cloud radius. This puts a 
lower limit on the cloud size according to Eq (1). If we assume that the photodissociated 
gas, as traced through its [CII] and [01] line emission forms a surface layer on the molecular 
cloud interiors, then we can estimate the cloud sizes. As discussed in Section 4.1.1, such a 
layer would have a thickness such that Ay ~ 3, and the atomic to molecular gas mass ratio 
requires that the atomic surface constitute ~ 3 x 10~ 2 of the total cloud mass. For spherical 
clouds, this implies a cloud radius corresponding to Ay ~ 300. Converting to physical sizes, 
these arguments then give T c \ ou d ~ 2.0 pc, at our modeled density and solar dust abundance. 
In a plane-parallel geometry, appropriate for clouds illuminated from a single direction, the 
geometric argument would suggest a cloud thickness corresponding to A v ~ 100, or 0.7 pc. 
Adopting the 2 pc spherical clouds for argument, each has a mass of 1.1 x 10 5 M , and there 
about 330 such clouds in the central 180 pc region. For reference, we note that increasing 
the CO abundance does not change the arguments above regarding the individual cloud 
geometry, but would correspond to a smaller molecular gas mass, lower volume filling factor, 
and fewer number of clouds. Table 2 summarizes the parameters for our adopted model. 

5. Conclusions 

We have detected bright CO J = 7 — > 6 emission from the nucleus of NGC 253. We 
combine our observations with lower-J measurements of both 12 CO and 13 CO and an LVG 
model to estimate the CO column density and physical conditions of the molecular gas in the 
starburst nucleus. The column density of CO we infer based on the transitions above J = 1 — > 
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is slightly larger than but similar to the estimates made by HHR99 for the central 23". 
However, we find that the excitation conditions of this gas must be very high to reproduce 
the observed J = 6 — > 5 and J = 7 — > 6 intensities. A single component with T ~ 120 K, and 
n ~ 4.5 x 10 4 cm~ 3 provides the best fit to all the observed transitions above J = 1^0 for 
both 12 CO and 13 CO, and is consistent with the ISO molecular H2 observations. We adopt 
this as the favored model. While a smaller additional mass of cool gas is required to account 
for the 13 CO J = 1^0 intensity, and more cool gas may be present, the CO column density 
of warm molecular gas cannot be much less than our derived values. 

The large mass of warm gas is very difficult to explain with a photodissociation region 
scenario. The molecular gas mass is a factor of 10-30 greater than the atomic gas mass, 
while the two should be comparable if heated with stellar UV as in the Galactic PDR 
templates. Dynamical heating of molecular gas from shocks is also unlikely based on energetic 
considerations, and given the lack of evidence for significant heating in fast shocks. On the 
other hand, the very large cosmic ray flux in the starburst nucleus of NGC 253 readily 
provides the necessary energy to heat the gas, and a means to deliver energy throughout the 
volume of large, UV-shielded clouds. It is likely the dominant heating term in the molecular 
ISM energy budget. While some of the warm molecular gas that we observe is likely heated 
by the far-UV radiation and shocks, our results strongly support the conjecture by (Suchkov, 
Allen, & Heckman 1993) that the bulk of molecular gas in starburst nuclei is held at high 
excitation by the enhanced cosmic ray flux. 
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Fig. 1. — CO J = 7^ 6 spectra toward the nucleus of NGC 253. The (0, Opposition is at the 
IR nucleus at, R.A = h 45 m 06 s .O, Dec. = -25° 33'40" (1950) and the beam size is 11". 
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Fig. 2. — LVG model calculations compared with observed line ratios. Contours of the 12 CO 
J = 7 ^ 6 to J = 3 ^ 2 intensity ratio are plotted as solid lines (values of 0.3, 1, 3, 10). 
Calculations are also performed for the 13 CO transitions assuming an abundance ratio of 50. 
The 12 CO to 13 CO J = 3^2 ratio is shown with the dashed lines (values of 5, 10, 20, 40). 
For both ratios, the observed values are shown with thick lines, bracketed with thin dotted 
lines showing 20% uncertainty. Scaling the CO abundance from the canonical 8 x 10~ 5 is 
equivalent to adjusting the velocity gradient in the opposite direction with the same factor. 
With the canonical CO abundance, the velocity gradient must be greater than 10km s -1 pc _1 
to reproduce the data. Larger CO abundances require even larger velocity gradients. 
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Temperature (K) 

Fig. 3. — Relationship between temperature, velocity gradient, and column density required 
to reproduce the observed intensities. At each point, the density is determined to match the 
12 CO J = 7^6 to J = 3^2 line ratio, then the observed integrated intensities are compared 
with the LVG model to determine the column density required. Higher temperatures and 
lower optical depths require smaller column densities. The constraints based on the 13 CO 
to 12 CO J = 3^ 2 ratio are shown as boxes for dv/dz= 10, 20, 40, and 80 kms _1 pc _1 . The 
sizes of the boxes are estimates of the uncertainties. 
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Fig. 4. — Integrated 12 CO and 13 CO line intensities in the nucleus of NGC 253 with our 
adopted model of the excitation. The data are compiled from the literature and smoothed 
or adjusted to correspond to a 15" beam, see Table 1. The solid lines are the total, the 
thin dotted lines in the lower left show a diffuse component which is energetically negligible 
except for a contribution to the J = 1 — > lines. 



